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We study the possibility to induce an effective gauge field for light confined to a Kagome lattice 
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that the Haldane Quantum Hall effect arises simply through a site-dependent phase (but constant 
amplitude) of the dynamic modulation. Within this scheme, we further demonstrate the existence of 
topological one-directional edge states immune to back-scattering losses, and discuss the possibilities 
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Topological order has opened a new frontier in the clas¬ 
sification of distinctive phases of matter, and is thus a 
center of attention of theoretical and condensed matter 
physics pQ. Its study has also reached the field of pho¬ 
tonics for two main reasons. First, photonic ana¬ 
logues of topological systems are a promising route to 
bridging theory and experiment. Second, a signature of 
a topologically non-trivial material is the presence of one- 
directional edge states providing energy transport im¬ 
mune to disorder. This could prove extremely valuable 
for slow-light photonic devices, which find a variety of ap¬ 
plications mi, but whose performance is severely lim¬ 
ited by back-scattering due to fabrication imperfections 

mi- 

Historically, topological order was first recognized in 
relation to the Quantum Hall effect. In that area, Hal¬ 
dane had a ground-breaking contribution in demonstrat¬ 
ing that the effect can arise even with zero magnetic 
field averaged over a primitive cell. The research into 
topological photonics was also started by Haldane in two 
theoretical studies m l, which were quickly followed 
by an experimental realization of a photonic topologi¬ 
cal insulator using gyromagnetic media [TO] , This re¬ 
sult was however obtained in the GHz frequency range. 
Due to the lack of suitable materials, reproducing this 
scheme in the visible or the near-infrared spectrum - 
which are the most interesting for applications - is still 
a major challenge. The milestone of an experimental re¬ 
alization of topological edge states for light in the near- 
infrared has been reached using coupled microring res¬ 
onators nn or coupled waveguides m by taking advan¬ 
tage of the symmetry-induced degeneracy of rotating and 
counter-rotating modes. More specifically, these systems 
are characterized by a preserved time-reversal symme¬ 
try (TRS), which leads to an important limitation of 
the topological protection. The ground-breaking result 
(which is now known as the Spin Quantum Hall effect) 
of Kane and Mele m that, for electrons, this protection 
is still present in TRS systems, relies on the anti-unitarity 
of the time-reversal operator (T 2 = —1). For photons, 



FIG. 1. (a): Kagome lattice with three sites in the primitive 
cell, and the corresponding Brillouin zone. Effective magnetic 
flux through a hexagon, similar to the one of the original 
Haldane model, is used here to open topological band gaps, 
(b): The corresponding band-structure with first-neighbor 
coupling J and zero flux <f>. There are six Dirac cones, and in 
addition a flat band. 


this operator is unitary, and the result no longer holds 
[2], at least not in its full strength. Instead, the pro¬ 
tection relies on the symmetry that prevents the mix¬ 
ing of propagating and counter-propagating modes in a 
waveguide, which in practice may be broken by disorder. 
This suggests the need for systems where TRS is broken 
01EHT9]. Recently, the possibility to use a fine-tuned 
dynamic modulation of a system to engineer a gauge field 
for photons has been shown both theoretically and 
experimentally [20J [21]. This scheme is employed here 
to induce a Haldane-like magnetic flux for photons on a 
lattice of optical resonators. 

The seminal work by Haldane [22 considered a hon¬ 
eycomb lattice (see the Supplemental material) with real 
first-neighbor and complex second-neighbor couplings. In 
the absence of the latter, the band structure of the lattice 
has six Dirac points, and no band gap. Haldane showed 
that the complex second-neighbor hopping terms, which 
result in zero average magnetic field over the unit cell, 
but non-zero magnetic flux through a triangle enclosed 
by second-neighbor hopping, break the TRS and open 
a topological band gap. Recently, this was successfully 
observed in a system of cold atoms in a ‘shaken’ opti¬ 
cal lattice [23] . which, together with previous research 
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FIG. 2. (a): Quasi-energy bands computed through diagonalization on the Floquet basis, for J — 0.1D, Ao = 0.9S2, ip = 2.1. 
The bands are repeated in orders of with m an integer, (b): Zoom-in on the m — 0 region of (a), (c): Bands computed 
through a perturbative expansion of the effective time-independent Hamiltonian. In (b) and (c), the Chern number for each 
band is indicated. 


in that field [24H26] , inspired the results presented here. 
In this work, we show how an analogue of the Haldane 
model can be achieved in a Kagome lattice of photonic 
resonators using a time-periodic modulation of the reso¬ 
nant frequencies, where only the phase of the modulation 
varies among different sites, in a spatially periodic man¬ 
ner. We further show the existence of back-scattering- 
immune edge states, and discuss the possibilities for a 
practical implementation of the system. 

Model- We consider a lattice of optical resonators, in 
which the resonant frequencies Ui are subject to a peri¬ 
odic modulation in time. The linear photonic Hamilto¬ 
nian, most generally, reads 

H = + Ai cos (fit + (f>i))a\ai — Jij a\ aj , (1) 

i ij 

where a t is the photon creation operator, J\j are the hop¬ 
ping coefficients, and Ai and fa denote the site-dependent 
amplitude and phase of the dynamic modulation, which 
can be achieved for example through electro-optic modu¬ 
lation Ea ns, optically-induced material non-linearities 

[29] , or optomechanical interaction with phonon modes 

[30] . The Hamiltonian is particle-number preserving, 
thus eq. 0 describes the system with any fixed num¬ 
ber of photons (sub-spaces of different photon numbers 
are decoupled). The equation also applies to classical 
light, since it is a concise way to write the coupled-mode 
theory that can be used for an array of optical resonators. 
In the Supplemental material, we outline the theoretical 
details of the Floquet theory [3TH33] that we employ to 
solve the time-periodic Hamiltonian of eq. 0. The way 
the Haldane model arises is best revealed through Flo¬ 
quet perturbation theory. There, to first order in 1/Q, 
new effective couplings can be derived: Jij -T J[- + iJ-j, 
with Jh , J"j real constants (see the Supplemental mate¬ 
rial). In our system, the imaginary iJ['■ thus introduces 
the magnetic flux required for the Haldane effect. 

The most straightforward way to achieve this effect 
would be by replicating the system of Ref. [23] through 


an appropriate modulation of a honeycomb lattice of res¬ 
onators. However, as discussed in the Supplemental ma¬ 
terial, this would require a spatial gradient of the modu¬ 
lation amplitude A{. This breaks the spatial periodicity 
and makes it impossible to analyze the system in momen¬ 
tum space, which is a significant theoretical disadvan¬ 
tage. In addition, in view of potential experimental real¬ 
izations, this feature introduces an extra challenge, since 
the maximum amplitude of the modulation is inevitably 
limited, which in turn would limit the maximum system 
size. Fortunately, this can be easily overcome through a 
modification of the lattice geometry - namely, by consid¬ 
ering the Kagome lattice illustrated in Fig. g a). This 
lattice has three lattice sites per elementary cell, and the 
band structure (Fig. 0b)) is similar to the one of the 
honeycomb lattice in that there are six Dirac cones. The 
main difference comes from the additional flat band. Im¬ 
portantly, in the presence of a flux similar to the one of 
the Haldane model, topologically non-trivial band gaps 
can be opened between the first and the second and/or 
the second and the third bands [34, 35.. 

The Kagome lattice studied here has identical res¬ 
onators of frequency cjo on all sites, and first-neighbor 
couplings only (along the black lines of Fig. 0a)), with 
a hopping coefficient J. The dynamic modulation that 
we assume has the form uja = c^o + Zlo cos(Qt + cp), 
ujb = + A) cos (Qt + 2(p), ujc = <^o + A) cos(Qt + 3(/?), 

where A, B, and C refer to the three sites of the primi¬ 
tive cell. Since the modulation is time-periodic, for time- 
scales larger than the period T = 27r/Q, we can apply 
the Floquet theory of quasi-energies [3TH33] , which is 
outlined in detail in the Supplemental material. There, 
we give two possible ways to compute the Floquet spec¬ 
trum of the system. This spectrum has a Brillouin-zone- 
like structure in the sense that, if uj is an eigenstate of 
the time-dependent Hamiltonian, so are uo + raQ, for all 
integer m. The eigenstates are time-periodic functions 
with period T, and can be expanded on the Floquet ba¬ 
sis \{ni},m ), where ni is the occupation number of site 
i. This results in an infinite-dimensional matrix for di- 
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agonalization with matrix elements given in the Supple¬ 
mental material. We consider the sub-space of a single 
excitation only, i.e. m = 1, rij^i = 0, and truncate the or¬ 
ders of m at ra max = 10 (convergence is always checked), 
which yields a finite matrix that can be diagonalized nu¬ 
merically. 

The Floquet two-dimensional band structure of the lat¬ 
tice can then be computed in this way, with time- and 
space-periodic solutions 

u n (k,t) = J2vi,m(Kn)e- ikR <e imQt , (2) 

i,m 

with Vi jTn ( k, n) the eigenvectors from the diagonalization, 
and Hi the position of site i. The Floquet band diagram 
is shown in Fig. [2ja)-(b) for J = O.lfl, A 0 = 0.911, 
p = 2.1. As discussed and displayed in panel (a), the 
bands are repeated in frequency space at an interval of Q. 
In panel (b), which shows a close-up of the zero-th order 
bands of panel (a), we see that band gaps are opened due 
to the dynamic modulation. To quantify their topological 
properties, we compute the Chern number for all bands 
by integrating the Berry curvature ^(k) |36l (37j over 
the Brillouin zone. Numerically, we compute ^(k) on a 
discrete mesh in k-space using the eigenvectors v m ^(k, n) 
[33135]. The non-zero Chern numbers (1, -2 and 1 for the 
three bands, respectively) confirm the non-trivial nature 
of the band gaps. 

The second way to handle eq. 0 (see the Supplemen¬ 
tal material) is through a perturbative expansion for an 
effective time-independent Hamiltonian. As mentioned 
above, this has the advantage of making the connection 
between this system and the Haldane model manifest, 
since the first-order terms in the expansion are imagi¬ 
nary couplings that introduce a flux in the red triangle 
of Fig. [ljb) . In Fig. |2|c), we show the bands computed 
by diagonalizing this effective Hamiltonian, which agree 
very well with the exact solution of panel (b), and the 
computed Chern numbers are the same. 

Edge states- Topological invariants like the Chern 
number cannot change as long as the band gap remains 
open. Hence, the width of the band gap is an impor¬ 
tant parameter, giving an energy scale to the topological 
protection against disorder (only fluctuations on a larger 
scale can destroy the topological properties). Thus, in 
Fig. [3j we plot maps of the gap width A t (if two gaps are 
present, the largest value is taken), versus the parameters 
Ao and ip. The data in panels (a) and (b) are computed 
for J = O.lfl, with the perturbation theory Hamiltonian 
in (a), and the full diagonalization in (b), and show very 
good agreement. In panels (c) and (d), J = 0.511 was 
used, and the agreement is no longer present. It is natu¬ 
ral that the perturbative expansion works well for small 
J/fl when the Floquet bands of different orders are well- 
separated (Fig[2ja)), but has limited reliability as J in¬ 
creases. Importantly, however, the topological effect is 
present even beyond perturbation theory: a gap of width 


Ay/fi A t/Q 



0 2 0 2 

<p [rad] ip [rad] 

FIG. 3. The (largest) width of the opened band gap due to the 
dynamic modulation of frequency fl vs. the amplitude Ao and 
the phase angle p for the Kagome lattice with first-neighbor 
coupling J — O.lfl, (a): Floquet perturbation theory; (b): 
expansion on the Floquet basis, (c)-(d): Same as (a)-(b), but 
for J — 0.5fl. The color scheme is the same in panels (a) and 
(b), as well as in panels (c) and (d). 


larger than 0.211 is opened for J = 0.511, Aq = 1.611, 
ip = 2.1. Notice that for any value of the parameters 
in this system, the band gap is inevitably limited to a 
fraction of II due to the higher-order Floquet bands. 

In Fig. [4] we show the band structures with the largest 
band gaps for J = 0.3fl, J = 0.5fl and J = 0.7fl, with 
parameters Aq and p chosen for the largest A t (see Fig. 
2 of the Supplemental material). Topologically, there is a 
difference between the bands in Fig. [2|b) and[4ja), with 
Chern numbers 1,-2, and 1, and those of Fig. 0b)-(c) 
with Chern numbers 1, 0, and -1. What is important, 
however, is that in both cases there are bands with a 
non-zero topological invariant. The bulk-boundary corre¬ 
spondence principle BE] then applies, guaranteeing the 
existence of gapless edge states at an interface between 
the topological material and a topologically trivial one 
(e.g. empty space). In terms of practical applications, 
propagating modes robust to disorder are thus expected 
to appear in a finite system. 

The existence of the topological edge modes is illus¬ 
trated in Fig. [5] for a ribbon geometry, with a finite 
number of sites in one direction, and periodic boundary 
conditions in the other. The one-dimensional Floquet 
band structure can again be computed by expanding on 
the Floquet basis, and is shown in panel (a) and (d) for 
J = 0.511, Ao = 1.611, p = 2.1. The difference between 
the two panels comes from the truncation at the edges - 
compare panels (b) and (e). Regardless of how we trun¬ 
cate, there is a band that closes the band gap of the 
bulk structure, due to the non-zero topological invari- 
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FIG. 4. The bands structure with the largest possible band gap for various values of J/Q. (a): J — 0.3Q, Aq — 0.50, ip — 2.1; 
(b): J — 0.50, Aq — 1.60, (p = 2.1; (c): J — 0.70, Aq = 3.050, p — 2.67. The Chern number for each band is indicated. 


ants. Modes belonging to that band are localized close 
to the boundaries of the ribbon; the important point, 
however, is that the modes at k x and — k x are localized 
at opposite edges. This is illustrated in panels (c) and (f), 
where we plot the position dependence of the magnitude 
of the eigenvectors of the two states indicated by a blue 
and a red dot in panels (a) and (d), respectively. The 
amplitude on the x-axis is the quantity ^ |r> m ^(k,n)| 2 , 
where the sum is over all m, and over all sites at the same 
position along y. The edge modes are exponentially lo¬ 
calized at the boundaries (notice the logarithmic scale 
on the x-axes of panels (c) and (f)), thus the overlap 
between the forward and backward-propagating modes 
decreases exponentially with the width of the ribbon in 
the y-direction. This is only possible due to the broken 
TRS, and ensures protection against back-scattering in 
the presence of disorder. 

Discussion- Several considerations have to be made 
for the results presented here to have practical implica¬ 
tions. We have not considered the loss rate ft of the 
optical resonators, which is in practice always non-zero. 
To be able to meaningfully talk about light transport, 
this must be smaller than the coupling constant J. In 
addition, ft must also be smaller than the band gap At, 
so that the latter can be resolved. By extension, this 
also implies ft <C In state-of-the-art photonic crys¬ 
tal cavities, k/ljo of the order of 10 -6 can now be rou¬ 
tinely achieved [4QH42] at telecommunication frequencies 
ujo/2tt « 200THz, thus k/2tt = 0.5GHz is a reasonable 
and conservative assumption. The coupling constant J is 
the easiest parameter to control by varying the distance 
between resonators. Thus the more important challenge 
is to have a sufficiently high A T . In fact, independently 
of ft, At is a general figure of merit for the magnitude of 
the topological protection that should be maximized. 

In Ref. [IT], electro-optic modulation was suggested as 
the practical tool for driving the resonant-frequency os¬ 
cillation. This offers sufficient control over the phase, and 
has been shown to be scalable Hung. The maximum 
achievable modulation frequency D/2i r is of the order of 
several GHz. A band-gap At of the order of 1GHz could 
thus be achieved, which lies just above the limit set by 




FIG. 5. (a): Floquet bands for the ribbon geometry shown in 

(b) , with a finite number of sites in one direction (the system 
is truncated at the solid black lines), and periodic boundary 
conditions in the other (along the dashed black lines). The 
parameters are as in Fig. f3^b) : J = 0.5S4, Aq — 1.6f2, <j) — 2.1. 

(c) : The spatial dependence of the eigenstates marked in blue 
and red, respectively, in panel (a). The y-axis is aligned with 
the y-axis of panel (b). (d)-(f): Same as (a)-(c), but for a 
different truncation (compare (b) and (e)). 


ft. We note that this challenge holds both for our pro¬ 
posal and for that of Ref. m- Very recently fill], it was 
suggested to use the coupling of the optical resonators 
to localized phonon modes to induce the frequency mod¬ 
ulation. In this scheme, Q is fixed by the phonon reso¬ 
nant frequency, which can be as high as D/2i r = 10GHz 
in two-dimensional optomechanical crystals [44]. This is 
sufficiently large for our scheme, and the required phase 
control can be easily implemented through the phase of 
the lasers driving the mechanical oscillations mm - We 
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note that, when compared to Refs. m and m , our pro¬ 
posal has a significant structural advantage, as it involves 
identical resonators with no intermediate (link) elements. 
Another recent optomechanical scheme m investigated 
the Kagome lattice of resonators, focusing on creating 
and probing topological states for sound (i.e. phonons). 
Within that proposal, it is also possible to create topo¬ 
logical states of light, but the size of the band gap is 
shown to be proportional to the phonon hopping coeffi¬ 
cient. This is typically orders of magnitude smaller than 
the phonon resonant frequency, and thus also smaller 
than the best optical loss rate k that could possibly be 
achieved in state-of-the-art photonic devices. 

While both of the modulation schemes mentioned 
above could be employed for an experimental realization 
of our system, a third option is also worth mentioning. 
Using the optically-induced Kerr nonlinearity, repeated 
switching at a THz rate has been recently demonstrated 
in a micropillar cavity [29]. The maximum amplitude in 
such a scheme is limited to only a fraction of U, but as¬ 
suming 0/27r = ITHz, Aq = 0.050 (which can be read 
out of the sine-like dependence of the cavity resonant fre¬ 
quency measured in Ref. ESI). J - 0.20, and ip = 2.1, 
we obtain for our Kagome lattice a topological band gap 
of width 0.0330, i.e. At/27t = 33GHz. This value is very 
similar to the magnitude of the disorder-induced fluctu¬ 
ations in the resonant frequencies of nominally identi¬ 
cal photonic crystal cavities nang, and, furthermore, 
the latter can in principle be reduced by post-processing 
techniques ga Eg. The predicted A t is thus, on one 
hand, two orders of magnitude larger than the loss rate 
of state-of-the-art cavities, and on the other high enough 
to ensure a truly sizable protection against disorder. 

Conclusion- In conclusion, we have described a 
straightforward implementation of the Haldane-like 
Quantum Hall effect for light in a lattice of optical cav¬ 
ities, with an effective gauge field produced through a 
time-periodic modulation of the resonant frequencies. 
The site-dependence of the phase of the modulation 
breaks time-reversal symmetry and opens topologically 
non-trivial band gaps, which, in a finite geometry, yields 
propagating, back-scattering-free edge states. These can 
find applications for high bit-rate storage [4], for en¬ 
hanced non-linear effects e.g. for frequency conversion 
or generation of non-classical light for quantum informa¬ 
tion processing [49l [52] . and for enhanced radiative cou¬ 
pling between distant quantum dots for on-chip quantum 
computation [53 [54]. 

This work was supported by the Swiss National Sci¬ 
ence Foundation through Project N- 200020.149537. We 
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